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Abstract 

We generalize the technique of linked cluster expansions on hypercubic 
lattices to actions that couple fields at lattice sites which are not nearest 
neighbours. We show that in this case the graphical expansion can be ar- 
ranged in such a way that the classes of graphs to be considered are identical 
to those of the pure nearest neighbour interaction. The only change then con- 
cerns the computation of lattice imbedding numbers. All the complications 
that arise can be reduced to a generalization of the notion of free random 
walks, including hopping beyond nearest neighbour. Explicit expressions for 
combinatorical numbers of the latter are given. We show that under some 
general conditions the linked cluster expansion series have a non-vanishing 
radius of convergence. 
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1 Introduction 



Linked cluster expansions are series expansions about completely disordered lattice 
systems Jl|, p[. Under very general conditions on the action the high temperature 
series can be shown to have a no n- vanishing radius of convergence 0. In many 
cases the singularity closest to the origin can be related to the physical singularity 
of the free energy, that is, the phase transition. To achieve valuable information 
from the series representations one needs to compute high orders (for most recent 
reference cf. 0, |5|, |^, 0]). This is best done by means of graph theoretical methods. 
Usually the expansion is restricted to contact terms that couple fields located at 
nearest neighbour lattice sites. On the hypercubic lattice, which we refer to in the 
following, this leads to a considerable reduction of the number of graphs. The most 
striking property in this respect is that they have to have an even number of lines 
in every loop. 

Here we discuss the question of how to take care of interactions beyond nearest 
neighbours in an optimal way such that the classes of graphs have to be modified 
minimally. This implies to keep them reliably small in number. As will be shown 
below, the graphs of nearest neighbour interactions are sufficient for the general case 
also. The only complication that arises can be absorbed completely by a general- 
ization of random walks, for which explicit expressions are available. 

Let A denote a D-dimensional hypercubic lattice, either infinite A = IP or finite 
A = x^o^Z/LjZ), with all L i even (for convenience) and with periodic boundary 
conditons imposed. To be specific, here we discuss models that are described by the 
partition function 

N 

Z(J,v)= / Hd N $(x) exp(-5($,t;) + ^^J a (x)<l> Q (x)), (1) 

xeA icgA a=l 

where $ denotes a real, iV-component scalar field, and J are external sources. The 
action is assumed to be of the form 



1 N 

S($,v) = Y,V(*(x)) + ~ ^^>ab(x,y)^ b (y), (2) 

x x=£y£A a, 6=1 

with the lattice site action V($), supposed to be 0(N) invariant and appropriately 
bounded from below. 

Truncated correlation functions are obtained by differentiation of the generating 
functional 

W(J,v) = \nZ(J,v), 
W^ a2n ( Xl ,...,x 2n \v) = <$ ai (x 1 )---$ a2 „(x 2 „)> c (3) 

Q2n 



dJ ai (xi) • ■■dJ a2n (x 2r 



-W(J,v) 



J=0 
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Susceptibilities are denned as zero momentum correlations such as 



KbX2 = < $ a(x)$ 6 (0) > C (4) 

x 

D-1 

S a , b fit = ^T(]T xl) < $ a (a:)$ fe (0) > c . (5) 



X 4 = 



The linked cluster expansion is the Taylor expansion of the generating functional 
W(J,v) with respect to v(x,y), 



W(J,v)= (exp^^^,2/)^-^yKW^ 



(6) 

v=0 



It generalized in the obvious way to connected correlation functions according to 
(H). Multiple derivatives of W with respect to v(x,y) are managed by the identity 

n 

dW 1 / d 2 W dW dW \ 

+ "nTTTT "oTTTa ■ ( 7 ) 



dv ab (x,y) 2 \dJ a (x)dJ b (y) dJ a (x)dJ b (y) i 

In most cases under consideration, the hopping term of the action is assumed to 
couple fields only at nearest neighbour lattice sites, that is, 

, x ( 2k 5 ab , x = y ± u for some u, /oN 
V ^ = { '0 i otherwise ^ 

/x denotes the unit vector in the /ith direction, fj, = 0, . . . , D — 1. Correlation func- 
tions and susceptibilities become series in the hopping parameter k, with coefficients 
depending on the toplogy of the lattice and the interaction V. If V satisfies appro- 
priate bounds, the series can be shown to have a non-zero radius of convergence, 
and convergence is uniform in volume 0. 

In order to manage the complexities which increase rapidly with increasing order, 
one is led in a natural way to a graphical device. Truncated correlation functions 
get a representation as a sum over connected graphs, each one endowed with the 
appropriate weight. In this description a line of a graph represents a nearest neigh- 
bour contact term. It implies the constraint that the two attached vertices have 
to be placed at nearest neighbour lattice sites. In this way, selflines are excluded, 
but otherwise no additional exclusion constraints hold. This leads to considerable 
simplifications and fast methods to compute weights of graphs. In particular, on 
the hypercubic lattice, for the weight of a graph to be non-vanishing, every occuring 
loop in this graph is subject to the constraint that it has an even number of lines. 
In turn, the number of graphs can be reduced in a considerable way. A graph with 
L lines contributes to the order k l . 

The simplications just mentioned do no more apply for more general contact 
terms, that is, if v couples fields at lattice sites with larger separations. A line 
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still represents an interaction v but otherwise no nearest neighbour constraint holds 
anymore. At the first sight this amounts to a considerable enlargement of the class 
of graphs to be considered. In turn all the algorithms developed for the construction 
of graphs to high orders have to be extended as well. 

In this letter we show that the arising complications do not have any influence 
on the graphical expansions. Hopping parameter expansions for lattice models with 
pair interactions over larger distances can be cast into graph classes identical to 
those already known from pure nearest neighbour interactions. This implies that 
the same algorithms to generate all contributing graphs apply to this case. The only 
complication that arises is that the notion of free random walks assigned to chains 
of vertices with two neighbours have to be generalized appropriately. 

In the second part of this paper we will generalize the proof given in namely 
that the radius of convergence is non-vanishing for linked cluster expansions of mod- 
els with nearest neighbour interactions, to the present case. It will be shown that 
under some general conditions the connected Green functions and susceptibilities are 
analytic functions in the hopping parameters for a large class of pair interactions. 

2 Graphical expansion. 

2.1 Some basic graph theory 

A graph is a sequence of two objects and two maps 

T = (£ r ,B r ,£ r ,$r), 

with Cr and Br 7^ disjoint sets, the internal lines and vertices of T, respectively. 
Er is a map 

E T :B T - {0,1,2,...}, 

v -> E r (v), (9) 

that assigns to every vertex v the number of external lines Er(v) attached to it. 
The number of external lines of T is given by Er = XXeB Er(v). Finally, $r is the 
incidence relation that assigns internal lines to their endpoint vertices. We consider 
lines as being unoriented, so $r maps onto unoriented pairs of vertices 

$r : C T -> (Br x B T ). (10) 

In order to stay more general we do not exclude selflines in this definition. External 
vertices are those that have external lines attached, 

B r , ext = {veB r \ E r (v) ^ 0}, (11) 

whereas internal vertices don't, Br,int = Br \ Br,ext- For every pair of vertices 
v, w G Br, m(v, w) denotes the number of common lines of v and w, i.e. the number 
of elements of 

M(v,w) = {I e C r | $ r (0 = (v,w)}, (12) 
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m(v,w) = \JA(v,w)\. 

w is called a neighbour of v if m(v, w) 7^ 0, and 

M(v) = { w G B r | m{v, w) ^ 0} (13) 

is the set of neighbours of v, with N(v) = \Af(v)\ its number of elements. For this 
definition, multiple lines are not count according to their multiplicity. For every 
integer n, n-vertices are those that have precisely n neighbours, 

B r n) = {v e B r | N(v) = n}. (14) 

An n-vertex has at least n internal lines attached. The number of total lines, both 
internal and external ones, attached to a vertex v is denoted by 

l(v) := m(v, w) + Er(v). 

weB r 

If l(v) is even for every v G Br, T is called even. In particular, then, T has an even 
number of external lines. A vertex v G B r is called weak n-vertex, v G i3p"^ eak , if 
N(v) = l(v) = n. A weak vertex does not have any external line attached. 

The following topological notions will be used below. A path is an ordered 
non-empty sequence p = (wi, . . . ,w n ) of vertices with m(wi,Wi+i) ^ for all % = 
1, . . . , n — 1. pis called a path in Y from w\ to w n of length n — 1. If in addition 
i« n = Wi, p is called a loop, with n — 1 its number of lines. If for every pair of vertices 
v,w G Br, v ^ w, there is a path from v to w, the graph T is called connected. A 
path p = (w 1 ,w 2 , ■ ■ ■ ,w n ) is called a maximal weak 2-chain if Wi G #p 2 ^ eak for all 

i = 2, . . . , n - 1, but w u w n & i3 r 2 ^ eak . 




Figure 1: Example of a maximal weak 2-chain of length 4. The 2-chain has 3 
vertices belonging to £>p 2 ^ cak and hence 4 lines. 

Two graphs 

T 1 = (d, B u E ± , $0, T 2 = (£ 2 , B 2 , E 2 , $ 2 ) (15) 

are called (topologically) equivalent if there are two maps fa : B\ — > B 2 and 2 : 
L\ — > £ 2 , such that 

$ 2 O0 2 = ^0$!, 

^0 0! = E U (16) 
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where o means decomposition of maps, and 



( f )l :B 1 xB 1 -> B 2 x B 2 

fa(v,w) = (Mv),Mw))- (17) 

A symmetry of a graph T = (£, £>, -E 1 , $) is a pair of maps <pi : £> — > £> and 2 : £ — » 
C, such that 

$ O 02 = 01 O 

= £. (18) 

The number of those maps is called the symmetry number of T. 

The set of equivalence classes of connected graphs with E external and L internal 
lines is henceforth denoted by Q E (L), and 



U Ge{L). (19) 



L>0 

If we impose the additional constraint that every loop has to have an even number 
of lines, we get the sets Q e ^{L) and . 

In practice one introduces further classes such as one-particle (or better to say: 
one-line) irreducible ones. Because this is of no importance for the following, we 
do not introduce them here. Also, just for simplicity we discuss the one-component 
model only. The generalizations are straightforward. 



2.2 Weights 

The graphical representation of the linked cluster expansion of susceptibilities is an 
expansion in term of equivalence classes of connected graphs as defined above. For 
every such equivalence class we need exactly one representative. 

Every graph T represents a number, which is called its weight. If T is connected 
and has E external lines, it contributes this weight to the susceptibility xe- For 
the one-component field models, it is computed along the following lines. Only even 
graphs need to be considered, otherwise the weight vanishes. 

1. All the vertices of T are placed at lattice sites. No exclusion principle holds, 

i.e. any number of vertices can be placed at the same site. By assumption, 
the model lives on the lattice A, which satisfies translation invariance. An 
arbitrary selected vertex is located at a fixed lattice site, avoiding a volume 
factor. A priori, all the other vertices of the graph are placed arbitrarily. 

2. Every internal line of T is assigned a factor v(x, y), where x, y G A are the lattice 

sites its endpoints are placed at. 



6 



3. The sum is taken over all possible placements of the vertices at lattice sites 

(except for the fixed vertex). The amount of computational work can be 
reduced considerably by taking into account the particular form of v(x, y) 
from the very beginning. 

° c 

4. A factor v 2n for every vertex v G Br, where 



v 



2n~ 



W {2n) {0,...,0\v = 0), (20) 



cf. Eqn. (^J), and 2n is the sum of internal and external lines attached to v. 

5. Two final factors complete the weight. The first one is given by (S'(r)) _1 , where 
S(T) is the symmetry number of T. It has the product representation 

S(T) = S P (T)- J] m(v,w)l, (21) 

(v,w)eBr xBr 

where the integer number Sp(T) corresponds to the symmetries of T under 
permutations of vertices. The latter product is due to the exchange of multiple 
lines between vertices. The other final factor counts the various enumerations 
of the external lines, 

E\ 

TW^RI- (22) 

Weighted susceptibilities such as the moment /12 are described by the same di- 
agrams as unweighted ones. The only difference is that they have the appropriate 
weight factors inserted. This, of course, must be done before the sum over the lattice 
imbeddings is carried out. 



3 Reorganization of graphs. 

Let us repeat that A denotes the D-dimensional hypercubic lattice A = x^^Z/LjZ, 
with all Li even, possibly infinity, and with periodic boundary conditions imposed. 
In terms of the classes of graphs introduced in the last section, truncated correlation 
functions and in particular susceptibilities allow for a representation as sum over 
graphs 

X e = Yl ^( r )' ( 23 ) 

with corresponding weight w(T). At this point we are not concerned with questions 
of convergence. Our interest here is in combinatorical rearrangements to arbitrary 
large but finite order. If desired we limit the discussion to propagators with v(x, y) = 
if I \x — y\\ > R for convenient R > 0, where 

D-l 

Il x ~2/H = Yl in |l x i _ Vi + nL t\, (24) 
i=0 
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and to graphs r G Qe with a limited number of lines. The limitation is not necessary 
for a convergence proof of the linked cluster expansion. This will be discussed in the 
last section of this work. We also mention in passing that the discussion below works 
also for more restricted classes of graphs such as 1-vertex irreducible or 1-particle 
irreducible ones. 

In the following we restrict attention to pair interactions that respect the eu- 
clidean lattice symmetries. A convenient notion in this respect is given by 

Definition 3.1 A signature p of length I G N is a collection of D non-negative 
integers l,eNU{0}, 

p= (Jo,--., to-i), (25) 

satisfying J^fjo 1 k = l- 

Collections of integers are unordered sequences, i.e. permutations of elements in 
p represent the same signature. For a signature p we write l p for its length, and the 
set of ordered D-tupels originating from p is denoted by Vd{p)- 

We may represent a signature p = (/o, • • • , Id-i) as a path connecting two sites 
x, y on the D-dimensional lattice. The path is given by {(xxi), (£1X2), ■ ■ ■ , (xo-iy)}, 
x i+ i = Xi + keii, i G {0, . . . ,D — 1}, &i G {—1, 1}, x Q = x, xp = y, where i is the 
unit vector in i-direction. 

We consider hopping propagators of the form 

v(x,y) = ^22K p v p (x,y), (26) 
pes 

where S is a finite set of signatures p of lengths l p > 1, and with 

leV D (p) e ,...,e £) _ 1 =±l 

5 T denotes the Kronecker delta on the torus, that is, 

8^ y = 1 if and only if inf \xj — yj + nLj\ = for all j = 0, . . . , D — 1, 

and 5^ y = in all other cases. Eqn. (p7|) denotes a propagator of a particular 
geometric profile. We assume that every p G S of the form satisfies for all 

./ 1 /; 1 

lj < mm — . 

J i=Q,...,D-l 2 

In particular, then, we have for p G S of length l p 

v p (x,y) = if \\x-y\\ +l p . (28) 

The assumption made here that S is a finite set is just for convenience. Sufficient 
convergence conditions that allow S to be infinite are discussed in the next section. 
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Our aim is to cast ( P3|) in such a form that we have to sum over graphs only that 
we have met already for the case of pure nearest neighbour interactions. That is, 

X e = E wCV ( r )> ( 2Q ) 

with appropriate weight w ev (r). Toward this end, let Qe denote the set of pairwise 
inequivalent graphs with E external lines and no weak 2-vertices, that is, 

Ge = {re G E \B?l eak = 0} . (30) 

We obtain from (|23l) 

Xe = E ^( r ), (31) 

where the weight is computed as outlined in the last section, the only modification 
being that a line now represents, instead of v(z,w), a factor 

oC r — ^ r — ^ qC - 

Pz-*w{{Kp)p<=S,V2) = 2^ 2^ v(z,Xi)v(x 1 ,X 2 ) • • ■v(x m - 1 ,W) (v 2 ) m ~ 

m>l xi,...,x m _i6A 

= ^(l-v 2 v)~\z,w). (32) 

The next step is to decompose P z -> w in terms of random walks of given length in 
units of lattice links. This is formulated as 



Lemma 3.1 For the hopping propagator v given by ([Mj), we have 

P z -,w((^p)peS, v 2 ) = P L w (( K p)peSi v 2 ), (33) 



with 

1 1 



^w,5 = E(^) m ~ E E 



m=l fi,-iPm6fi xi,...,a; m _i6A 



■v Pl (z, x x ) ■ ■ ■ v Pm {x m ^,w) (2KpJ • • • (2/t Pm ; 

' c\ E„ e s "p -1 



(34) 



E (S) ~ Msa^iKw. (35) 

("pipes) pG5 

where we have set xq = z, and with non-negative integers 



(36) 

K„ = 



Furthermore 



P z ^w(( k p)p&Si v 2) 7^ om y if I — \ \z — w\\ > Q and even. (37) 



Proof: Inserting ( Eg ) into the right hand side of the first equality of fl3"2| ) and 
collecting the sum over all products of propagators v p of fixed total length I = J2i ^p»j 
we obtain (|3~3"|), with Pj^^ as given by ([34]) . Furthermore, fl3~5f ) follows immediately 
from the Taylor expansion of 



-i 



Pz-*w((K p ) P es,v 2 ) = lc 1 1- v 2 ^2(2k p )v p \ (z,w). 

v 2 V pes J 

Finally, for every term of (Bl) we have, with xq = z and x m = w 



m—1 m—1 m—1 

\z — w\\ = || 

i=l i=l i=l 



and we have used ( p8|) and the triangle inequality for the norm (p4|). For j 
0,...,P>-1, 



m—1 



inf \zj — Wj + raL,-| = inf | > (x% — ^j-i), + nLA 

j=l 

m-1 , v 

|(xi — Xj_i)j + nLj\ J + even number. 



Hence, 



D-l 



\z — w\ 



Einf \zj — Wa + nLA 



j=0 
m—1 

= ||xj — Xi-i\ \ + even number 
i=i 

= I + even number, 

that is, (p7|). This completes the proof. □ 

Inserting the decomposition (p3|) into (|3TD , we obtain the desired representation 
(T2D|). The property ( |T7[ ) ensures that only those graphs contribute which have an 
even number of lines in each of their loops. 

Hence we end up with the same class of graphs as for the case of pure nearest 
neighbour interactions. The weight of a graph is computed as described in the last 
section with the following exception only. Every maximal weak 2-chain of a graph 
contributes a factor P z ^ w ((n p ), v 2 ) to the lattice imbedding number, where I is the 
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length of the chain and w, z the lattice sites where its two endpoint vertices are 
placed at. We notice that for even / neither these endpoint vertices have to be 
different nor the lattice sites z and w. If I — 1, P l z ^ w just imposes the nearest 
neighbour constraint on z and w. 

If we choose the hopping parameters (k p \p G S) according to 

for all signatures p G S, with the c p being fixed, every graph T G Q e ^ with L 
lines that contributes to ( |2§|) has a weight proportional to k l . In this case, every 
graph contributes to the same order of k as it would for the pure nearest neighbour 
interaction. 

Translation invariance amounts to Fourier transform. With 

n<£ £ ) (38) 

we have 

D-l 



eo,...,e£)_i=±l i=0 

The random walk numbers TVi^if 65 ^ are then given by 

(E„»,)! ^ 



U P np - A P e5 \/gP D ( P ) i=o / 

Similar as for the pure nearest neighbour interactions, this representation can be 
used to compute the random walk combinatorical integers Nhl peS) explicitly. 

As an illustrative example we supplement the nearest neighbour couplings by 
next-to-nearest neighbour interactions. In this case, we have 3 signatures, S = 
{Pi,P2,P3}, with 

Pl = (1,0,0) p 2 = (1,1,0) p 3 = (2,0,0) 



• • • 1 • o • 

written here for D = 3 dimensions. A maximal weak 2-chain of length I as part of a 
graph implies a factor Pz^ w ((Kp) p& s,v 2 ), cf. (|3^ ) ,(|35| ) , with z, w the lattices sites its 
end vertices are placed at. In order to get lattice imbedding numbers, one needs to 
compute the random walk numbers Mz^w 2 ^ for all nonnegative integers n\, n 2 , n 3 
with 

I = ni + 2(n 2 + n 3 ). 
In Table |l] we give all the numbers for z = w and chain length I = 16. 
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Table 1: Example of next-to- nearest neighbour combinatorical random walk num- 
bers Afz-^w 2 '™ 3 ^ on the 3-dimensional hypercubic lattice, infinitely extended in all 
directions. The numbers are for z = w and total chain length I = 16. The pure 
nearest neighbour case gives Af^Q^ = 27770358330. 



m 


n2 


"-3 




n\ 


™2 


"-3 


. /-(ni,n 2 ,n 3 ) 
■^0-0 


14 


1 





137889591840 


4 


4 


2 


20325816000 


14 





1 


56502916470 


4 


3 


3 


11722233600 


12 


2 





277007722992 


4 


2 


4 


3890527200 


12 


1 


1 


227358907776 


4 


1 


5 


678585600 


12 





2 


48295573326 


4 





6 


54299700 


10 


3 





288364684608 


2 


7 





437724000 


10 


2 


1 


355550656032 


2 


6 


1 


1266904800 


10 


1 


2 


151142397120 


2 


5 


2 


1617174720 


10 





3 


21791715660 


2 


4 


3 


1175005440 


8 


4 





165353348160 


2 


3 


4 


509120640 


8 


3 


1 


272212617600 


2 


1 


5 


142067520 


8 


2 


2 


173721477600 


2 


1 


6 


18506880 


8 


1 


3 


50101286400 


2 





7 


1610280 


8 





4 


5517544230 





8 





4038300 


6 


5 





51139740960 





7 


1 


13305600 


6 


4 


1 


105386339520 





6 


2 


20049120 


6 


3 


2 


89715669120 





5 


3 


17015040 


6 


2 


3 


38870314560 





4 


4 


9918720 


6 


1 


4 


8536207680 





3 


5 


2983680 


6 





5 


761454540 





2 


6 


1028160 


4 


6 





7704774000 





1 


7 





4 


5 


1 


19072972800 








8 


44730 



4 Convergence of the linked cluster expansion 

It was shown in B that the linked cluster expansion for connected Green functions 
and susceptibilities have a non-zero radius of convergence for the case of nearest 
neighbour couplings, i.e. v p = for all signatures with length l p > 1. We will 
generalize the proof given in || to hopping propagators which contain terms beyond 
nearest neighbour interactions. 

We do not suppose here that there are only a finite number of terms in the hop- 
ping propagator eq. (|26|) . The non-local hopping parameter may connect arbitrary 
lattice sites. S is the set of all signatures p of length greater or equal to one. 
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Let us define a hopping parameter k and constants c p for all signatures p by 



l p denotes the length of signature p. The constants c p have to be chosen such that 
the series on the right hand side of eq. (PB|) is convergent. The next lemma gives a 



sufficient condition on the coefficients c p . 

Lemma 4.1 Suppose that there exists a positive constant 7c > such that for all 
signatures p 

\c p \ < (2- 1 - r iK- 1 ) 1 " 7c, (40) 
where rj > is a constant obeying 

<AD2- D . (41) 

y.O^yeA 

Then a bound for the non-local part of the action is given by 

|~ £ $(x)v(x,y)$(y)\<4D\K\Y,mx)\ 2 (42) 

x,yeA xEA 

for all complex k, \k\ < 7c and all $ G TC(A) (=Hilbert space of square summable 
functions on A). 

Proof : The number of all signatures p with length l p = I is bounded by 

|{p signature] l p = l}\< 2 D+l . (43) 
Then the left hand side of inequality ( f42|) is bounded by 



E E mn^wi 2 . 



!p=l|a-l/ll 



Using supposition ( fl~0|) the bound becomes 

|/c| E E 2-*- 2/ " |$(x)| 2 . (44) 



By the bounds (fH) and (|44j) the left hand side of inequality (f42|) is estimated by 

2 D |/«| E 2-*-^" |$(a;)| 2 . 



a: ,y £ A 



Using supposition ([41]) we obtain the assertion. □ 

The convergence of the linked cluster expansion is presented in the following 
theorem for a symmetric translation invariant scalar field theory on an infinite lat- 
tice. A modification of the proof for iV-component models or models on a torus is 
straightforward. Also the restriction to translation invariance and Z2-symmetry is 
of no relevance for the convergence proof. 
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Theorem 4.1 Let a translation invariant and Z2- symmetric model be defined on 
the lattice A = TP by the partition function in eq. (flp. Suppose that the supposition 
of lemma is valid and there exist a positive c > and a real 5 such that for all 
$ £ R 

V($) > c$ 2 - 5. (45) 

Then there exists a positive constant such that for all x, y £ A £/ie 2-point Green 
function W™(x,y) is an analytic function in K(k*) := {k £ C| |ac| < ft*}. Further- 
more, there exist a constant a > and /or a// « £ t/iere zs m(«) > suc/i 
£/iat /or all x, y £ A 

|W (2 W)I <e a gf^ + aexp{-m||a;-y||}, (46) 



where 



(2) / d$$ 2 exp{-V($)} 



gy : Jd^$ exp{-y($)} ' ^ 
Furthermore, for k — > 0, 

= 0(ln |k|) > 0. (48) 

For m > the series in the definition @ of the 2-point susceptibility \2 is convergent. 
Thus, using a well-known theorem for analytic functions, the above theorem implies 
that the 2-point susceptibility X2 exists and is an analytic function for k in K{k^) := 
{k £ C\ \k\ < k*} for some > 0. 
Define a positive number 

q v ■= y"d$exp{-V($)} (49) 
and for all finite subsets X of A the partition function 

Z(X\J, k) = qy lXl I Y[ d$(x) exp{-S(X|$, k) + J ( X )H X )}, ( 50 ) 



where the action is 

5(X|$, K ) = ^y($(a;))-i *{x)v{x,y)$(y), 

xdX x,y£X 

and \X\ denotes the number of lattice sites of X. The following lemma states that 
the partition functions are well-defined if the local interaction V obeys the stability 
bound (E5|). 



Lemma 4.2 Suppose that there exist positive c > and real S such that the stability 
bound for the local interaction eq. ftffi ) is valid. Then we have for all k £ C and 
e > obeying 4d\n\ + e < c, and for all subsets X C A 

S(X|$, K )>e]T$ 2 (x)-5|X|, 
and the integral on the right hand side of Eq. fl5(\) is convergent. 
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This follows from the stability bound (23) and Lemma 3~T 



In contrast to the convergence proof of || given for nearest-neighbour couplings 
we need here a more general notion of polymers. Here we call all finite nonempty 
subsets X of A polymers. Define Mayer activities for polymers X by 

n 

M{X\J,K):=-5 im + Y,{-^Y~\n-l)\ II^I^)- ( 51 ) 

Then, the Mayer Montroll equation for the connected 2-point Green function reads 
(cf. H for a proof), for all x, y G A, 

w^(xy) = T T T d 2M(g|J;/t) 

Q: Q polymer n>0 Vi Vn£K(Q) V ' ViV 

(_ M )^ + -+^| J=0 (52) 
where we used the following definitions 

K(X) := {P = {Pi, • • • , P n }\ Pi C X, P a n P 6 = 0, Va ^ 6, n G N U {0}} 



and Po := {{ x }\ x G Q} i n the sum on the right hand side of eq. (|52|). Furthermore, 
for a set P = {Pi, . . . , P n } consisting of polymers P a we have defined 

ConniV) := {V G K(A)|VP' G V 3a G {1, . . . , n} : P' n P a ^ 0} 

and used the abbreviation 



,_ M) V 1+ - + V n ._ JJ(_ M (p.| JjK )). 



8=1 



The following lemma presents the conditions under which the Mayer Montroll ex- 
pansion eq. (|52|) is convergent and exponentially bounded (cf. H). For a proof see 



Lemma 4.3 Suppose that there exist positive constants a, «;* > and a function 
m(re) > stzc/i i/iai 

a 



£ |M(P| J = 0, «)| exp{a|P|} < - (53) 



and 



P:y&P 
\P\>2 



E r L/=o expaP < 

pf^P dJ{x)dJ{y) 
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e a <7y $Btf + ;r exp{-m||x -y||}, (54) 



for all x, y G A and \k\ < where 



(2) JtF** 2 exp{-^($)} 



9y : " /d"*exp{-V($)} • (55) 



T/ien the series expansion eq. ( pBj) is convergent and 

\W^(x,y)\ < e a q^6 xy + a exp{-m\\x - y\\}. (56) 

As n — > 0, 

m(K) = 0(ln \k\). 

The next lemma is the core of the convergence proof. It states sufficient conditions 
for the supposition of Lemma |4.3| to hold. 



Lemma 4.4 Suppose that the stability bound ($5) and the supposition of lemma 
are valid. There exist positive constants a > and > and a positive function 
m{K) > such that inequalities ( |5j[ ] and are valid for all \k\ < /«*. Furthermore, 
for k — > 0, 

m(«) = 0(ln 



Proof : For the proof we will use the tree graph formula [|1(J for Mayer activities. 
Let X be a polymer, \X\ > 2. Then 



M(X\J,k) 



dt X y\ 



11 (H>(x)v(x,y)$(y)) exp{- £ C™(*,2/) $0&M^2/)$(l/)}- (57) 

<xy>dT x,y£X, x^y 

In the following we will explain the notations used in formula fl57|). T(X) denotes 
the set of all tree graphs (graph containing no loops) with lines (xy), x G X and 
vertex set X. Furthermore, 

t™ n (x,y): = 

min {ti\ path connecting x and y and containing link / G r}. 



Using the stability bound ([45"D, definition ([49]) of the constant gy, Lemma [12] and 
the tree graph formula fl57|) we obtain 

\m(x\j = o,k)\ < E f II Hx,y)\)q v lxl e slxl 
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where 

d T (x) := \{l £ t\ line I emerges from vertex x}\ 

is the number of lines in the tree graph r which are connected to vertex x. For 
the estimation of the fields in the term n <xy>( z T Q(x)v(x,y)<&{i)) we have used the 
stability bound (flop and 



POD 

2 dx exp{— ex 2 }x n = e 2 
Jo 



■■ L ,n + 1. 



and we write z\ = r(^+l) for z > 0. For X = {x%, . . . , x n } and dj G N, i G {1, . . . , n} 
let T(X; di, . . . , d n ) be the set of all tree graphs r G T(X) such that d T (xi) = dt for 
alH G {1, . . . , n}. Cayley's Theorem ( for a proof cf. |fLl] ) is 



\W,-M\= (di _^...> dn _ 1)V (59) 



In eq. ( 58|) we replace the sum over tree graphs by 



E = E E 

r: 

Using that 



TireTpO di,...,d„>i reT(X;di,...,<Z„) 

E? =1 dj=2(n-l) 



(^)!<(d-l)! 

for all d > 3 and that there exists at least one iGl such that d T (x) = 1, we obtain 

u^^Y^y- ^ (^ ixhi n w*) - ( 6 °) 

Furthermore, 

\M(X\J = 0,k)\ exp{a|X|} 

|X|>2 

= (~ ~ S |M({y,x 2) ...,x„}|J = 0,/«)|exp{a;n}. (61) 

n:n>2 ^ «a.->«n 

We have replaced here the sum over all polymers X with \X\ > 2, y G X, by a 
sum over distinct labelled lattice sites X2, . . . , x n . We have defined xi := y. Then a 
tree r G . . . , x n }) corresponds to a tree T({1, . . . , n}) where the vertices are 

labelled by 1, . . . , n. We want to estimate the sum over all lattice sites x%, . . . , x n for 
a given tree r. For that we have to find a bound for 

j2 n w**>**)\ (62) 
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for all tree graphs r G T({1, . . . , n}). This can be done by estimating 

E \v(x,x')\. 



x'-.x'eA 



To obtain later in this proof the exponential factor exp{— m\\x — y\\} in the bound 
fl54] ) we find an upper bound for 



J\\f * \ \\x-x'\\ 

\k\ 



V |v(x,a;')|(|-7 

Using the suppositions (|40D and ( f41~|) and (|26|) , we see that for k with \k\ < < 7t 

E Hx,x')\(^-^ < E E 2 M«* p 



< 2 E E 2^2-*-* 

x':x^x'&k p:l p = \\x— x'\\ 

< 8Dk*. (63) 

We have used that Z p > 1 for all signatures p. We sum in expression ( |62"D successively 
over all X2, ■ ■ • , x n G A starting with vertices having a maximal distance to the "root" 
vertex y = x\. This gives a bound 

E II \< x ^)\ < (8D\K\) n -\ (64) 

Thus, using 

E 1 < 2 2(n - i} , 

tfl,...,Jn>l 
E|L 1 d i =2C»-l) 

Cayley's theorem Eq. (|59|), inequalities flBC|) and (|6T|) and the bound (|64"D, we obtain 
E |M(X|J = 0,«)| exp{a|X|} 

|J¥|>2 

^ E ^r(32(^)! V^^ex P {(a + 5)n} 

n: n>2 

= E -( 32 (^) ! D^qyh-h^re^qv 1 ^ = ~ ~ «) e a+ V^, (65) 

n>l 

where 

The proof of the bound of the term on the left hand side of ([54]) is similar to the 
proof of inequality (|65|). To obtain the exponential factor exp{— m\\x — y\\} we use 
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the fact that in the trees r occuring in tree graph formula of %^)dj(y) there exists 
a path connecting the vertices x and y. Then we use the bound 



n Hx a ,x b )\<cfy^ n 

(ab)er ' (ab)&T 



V(X a ,X b 



' K * \\\Xg-X b \\ 



For the sum over all positions X2, ■ ■ ■ ,x n for a given tree r we use inequality (|63l). 
The result is 



d 2 M(X\J,K) 



-_ \ exp{a\X\} 



X:y£X 
|X|>2 



where 



dJ(x)dJ(y) 
< — ln(l — u) e a+5 e~i exp{— m\\x — y\ 



m = 111(7-7). 

AC 



(66) 



This proves the assertion if ft* is small enough. □ 

Lemmata |4.3] and [4.4| prove that the Mayer Montroll equations ( |52D are abso- 
lutely convergent under the suppositions of theorem |4.1| . By a standard theorem for 
analytic functions we see that theorem [4.1| holds. 



5 Summary 

We have discussed the linked cluster expansion in application to lattice models with 
general pair interactions beyond nearest neighbour couplings. Under very general 
conditions, free energy, truncated correlation functions and susceptibilities are uni- 
formly convergent in volume to analytic functions, for sufficiently small hopping 
parameters. 

The graphical expansion is cast in terms of graph classes well known from the 
pure nearest neighbour interactions on hypercubic lattices. This implies that all the 
simplifications and computational efforts that are based on the topology of those 
lattices are fully exploited. The only modification occurs as a generalization of free 
random walks. 
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